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Abstract 

By employing retarded and advanced propagators, Aurenche and Becherrawy showed 
how to rewrite the real-time thermal Feynman rules so that the temperature de- 
pendence is removed from the free propagators and transferred to the vertices. The 
present paper introduces retarded and advanced field operators and incorporates all 
temperature dependence into the interaction term of the Hamiltonian. Perturbative 
expansions of the Green functions in the Hamiltonian formulation give the correct 
results order by order in perturbation theory. The spectrum of the temperature- 
dependent Hamiltonian contains the thermal quasiparticles that produce poles in 
the propagators. 
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Introduction 



The poles in the free, thermal propagator are at k 2 = m 2 . However in higher 
orders the poles occur at a temperature-dependent effective mass. This effect 
was first discovered in self-interacting scalar theories by Dolan and Jackiw [1] 
but is now known to be a quite general phenomena. The temperature depen- 
dence of the Higgs boson mass is responsible for the phase transition at large T 
in the electroweak theory. Temperature-dependent masses also occur for par- 
ticles with spin. In high-temperature QCD both quarks [2,3] and gluons [3,4] 
have effective thermal masses. The thermal masses are known to determine 
the location of the normal threshold branch points for particle production [5]. 

Although thermal mass effects are ubiquitous in perturbative computations 
of Green functions, it is not clear how to account for them in the operator 
theory since the Hamiltonian contains no temperature. As an example, for a 
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real scalar field x a typical Green function that is useful at finite-temperature 
is the full retarded propagator 



,-0E A 



D' R {x -y) = -i9(x - y ) £ -— (A| x (x), X (y) \A). 



Unlike the free retarded propagator, the full retarded propagator depends on 
T in a very complicated way. The Fourier transform, D' R (k), has poles and 
cuts in k whose locations are temperature-dependent. How these poles are 
related to the spectrum of the Hamiltonian is unclear. The explicit Fourier 
transform 

D W = + Pa - P B ) k !<tTL (^^) 

a? b k + E A - E B + ir) \ Zp J 



offers no indication that the singularities in ko are temperature-dependent 
since Ea and Eb are the zero-temperature energies of multiparticle states \A) 
and \B). 

The standard formulations of finite-temperature field theory do not help ex- 
plain how the singularities in k become temperature-dependent. The pertur- 
bative Feynman rules may be derived using either the path integral [6-8] or 
operator methods [9,10] and specify that the vertices are the same as at zero 
temperature and only the propagators are modified by T. The result for the 
real-time thermal propagator of a free scalar field is a 2 x 2 matrix 



A ab (k) 



(l + n)A R -nA A e° ko n(A R -A A ) 
e (0-*)*o n (A fl - A A ) nA R -(l + n)A A 



XI) 
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in which A R = l/(k 2 — m 2 + ik^e) and = A * R are the free retarded and 
advanced propagators, n = l/(exp(/3/c ) — 1) is the Bose-Einstein function, 
and a is a free parameter in the range < a < f3. The choice a = (3/2 
is common. Often (1.1) is expressed in terms of the Feynman propagator 
Ap = l/(& 2 — m 2 +ie) and its complex conjugate [6-8]. However the retarded 
and advanced expressions will be more useful in what follows. 

A novel formulation of the Feynman rules was developed by Aurenche and 
Becherrawy [11] and will be of central importance here. They showed that the 
Feynman rules could be rewritten so that the propagators are independent 
of T and the vertices contain all T dependence. This approach was developed 
further by van Eijck, Kobes, and van Weert [12,13]. The notation of [13] will be 
most convenient here. The basic idea it to express (1.1) in the skew-diagonal 
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form 



A A (k) 
K A R (k) 




V yb (k ). 



(1.2) 



The important feature is that the temperature dependence from the Bose- 
Einstein function is contained in the matrices V(ko). This diagonalization 
leads to a reformulation of the Feynman rules in which the free propagator 
is the center matrix in (1.2) containing Af>(k) and A^(A:). That propagator 
matrix is temperature- independent. The matrices V(ko) are incorporated into 
newly defined vertices which are temperature-dependent. The transformation 
automatically gives exactly the same momentum space integrals as the stan- 
dard Feynman rules since it is only a rewriting [11-13]. 

The purpose of this paper is to use (1.2) to construct an operator Hamilto- 
nian containing all the temperature dependence. The new Hamiltonian should 
produce exactly the same physical effects, and in particular the same n-point 
functions as those obtained by the usual methods of finite-temperature field 
theory. The Hamiltonian comes from a Lagrangian and the motivation for 
the choice of the latter comes from the path integral for finite-temperature 
Green functions. Sec 2 explains the choice of the Lagrangian density and the 
field variables <pn and <pA- The Lagrangian density contains time derivatives 
from first order through N'th order, where N is large but arbitrary. Sec 3 
defines the canonical momenta associated with the N time derivatives of the 
fields and constructs the full Hamiltonian operator. As preparation for what 
follows, Sec 4 quantizes the free theory with no interaction and no N'th or- 
der time derivative. Sec 5 quantizes the free theory including the N'th order 
time derivative and confirms that the interaction-picture propagators have the 
correct behavior in the limit N ^oo. Sec 6 applies the standard operator per- 
turbation theory using the evolution operator U(ti,t 2 ) to compute full Green 
functions as vacuum matrix elements of time-ordered operator products. In 
the limit N — > oo the perturbative expansion of the Green functions is exactly 
the same as the usual Feynman rules, thus proving that the temperature- 
dependent Hamiltonian describes the same physics. Sec 7 discusses the results 
and implications. 



2 Motivation for the Lagrangian 

This section will explain why the Lagrangian density is chosen to be (2.17) 
but will not prove that it is correct. The proof of (2.17) occurs in Sec 6, 
where it is shown that (2.17) leads to the correct perturbative expansion of 
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the Green functions. To motivate the Lagrangian it is useful to start with the 
path integral formulation of thermal field theory in real time [6-8]. In real-time 
formulations it is necessary to double the number of fields. Instead of one scalar 
field x ^ is necessary to use two, denoted x+ an d X— The finite-temperature 
generating functional is 

Z[j\ = J D X+ D X - exp[* J d 4 x(C c + j +X+ + j-X-)}, (2.1) 

where the contour Lagrangian is 

C c = \{d» X+ Wx + ) - \m\l - A( X+ ) 4 

- \{d,X-Wx-) + \W- + A(x-) 4 - (2.2) 

To absorb the ultraviolet divergences the parameters m and A should be the 
bare mass and bare coupling constant, which should then be written as the 
physical mass and coupling plus counter terms. As ultraviolet divergences play 
no role in what follows, these details will not be displayed. 



2.1 Choice of C 



oo 



It is straightforward to express the quartic interactions as functional deriva- 
tives of the free generating functional [6-8], Z 00 [j] = Det 1//2 (iA afe ) exp (il^fj]), 
where W[j] is determined by the free, thermal propagator A a i>: 



J d 4 xd A yja(x)A ab (x-y)j b (y) 
= ~ 1 iJ ^U-k)A ab (k)j b (k) (a,b = ±). (2.3) 

This is the usual starting point for deducing the Feynman rules. 

Matrix Transformations: The notation used will be that of van Eijck, Kobes, 
and van Weert [13]. The free thermal propagator matrix A ab shown in (1.1) 
has row and columns labeled by Latin letters a, b = ± as follows 



lab 



A ++ (k) A + 4k)" 
A_ + (k) A_4k) J 



(2.4) 
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(In other references [6-8] the rows and columns are instead labeled by the 
integers 1,2.) The center matrix in (1.2) will be written 

=(*»<*> *»W)J ° A ^»V (2.5) 
\A Aa (k) A AA (k) ) \A R (k) ) 

Thus two subscripts are used to specify the free retarded and advanced prop- 
agators A R (k) and A A (k) as matrix elements. The matrix relation (1.2) is 
written 

Aab(k) = Vj a {-h)A aP {k)V pb {h). (2.6) 

Latin letters a, b, c, . . . always denote ±; Greek letters a, /3, 7, . . . always denote 
R/A. The matrix V/3f,(k ) has rows labeled by R/A and columns labeled by ±. 
The explicit form of V^{k Q ) is contained in [11-13] and is discussed in Appendix 
A. As in [13] the inverse matrix V~ 1 (k ) will be denoted by U T (—k ) : 

U T (-k )V(k ) = 1. (2.7) 

The rows of Upb(ho) are labeled by (3 = R/A and columns by b = ± just the 
same as for Vpb(ko). These automatically satisfy 

U(k )a 3 U T (-k ) = a 1 V(k )a 3 V T (-k ) = a 1 . (2.8) 

Retarded and Advanced Sources: A new version of the generating functional 
arises if the propagator relation (2.6) is substituted into (2.3): 

w[j] = -\j ^j a (-fc)e a (-fco)A Q ,(fc)^(g Jfc (fc). 

This suggests that W be considered as a functional of two new sources Jr/a- 
The transformation from j± to Jr/a and the inverse transformation are 

Mk) = Vp b (k ) 3b (k) Jp{k)U Pc {-k ) = j c (k) (2.9) 
In terms of the retarded and advanced sources, W is 

W[J] = ~^\w^ J "(- k ) A °^Mk). 



V(k )U T (-k ) = 1 
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The two non- vanishing terms in the above integrand are equal since A^(fc) 
Aar(— k). The same result in coordinate space is 



W [J] =-J d 4 xd 4 yJ A (x)A R (x-y)J R (y). 



Retarded and Advanced Fields: The new sources Jr/a couple most naturally 
to a new combination of fields. To find this combination, it is only necessary 
to convert the coupling between the sources and the fields from coordinate 
space to momentum space and then use (2.9) : 



J d A xj c (x)xc(x) = J j c (-k)Xc(k) 
r d 4 k 

= J j^y M~ k ) U Pc(ko)Xc(k). 

The last line suggests that Up c (ko)x c {k) should be defined as a new field 0^, 
for (5 — R or A. The transformation from x± to 4>r/a and the inverse are 

Mk) = Up c (k ) Xc (k) Mk)Vpa(-ko) = Xa (k). (2.10) 



The source term is now expressed as 
f d 4 k r 

J J p(~ k ^p( k } = J d ^ x + Ja{x)4>r{x)}. 



Free Lagrangian: The Lagrangian appropriate for the 4>r/a is determined by 
examining the free contour action 

J C = / d 4 x [\(d, X+ ){d» X+ ) ~ \{d»X-){d»X-) - \m 2 X 2 + + \m\ 2 _\. 



In momentum space the action is 

J ° C ° = ^ / (^4 - m2 )Xa(-k) v 3 ab Xb(k). 

In terms of the new fields 4>r/a the free action is 

d 4 k 



J oo 



= l *lj2nY (^- m2 )^(-^)^(^o)^K(-^o)0 7 (A:). 



From the matrix identity (2.8) this simplifies to 
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J£> = § / ^(* 2 -™>*(-*)4A(*) 

In coordinate space the action may be written 7q = / d A x £00 where 

Ao = (dM(d^ A ) - m 2 <t> R 4> A . (2.11) 

The reason for the double zero subscript will be explained in Sec 4. The 
check ( " ) appearing in £ 00 is a reminder that it is the difference of two La- 
grangian densities. When quantized it will be a Liouville operator [14,15] that 
determines the time evolution of quantum amplitudes and of their complex 
conjugates (as suggested by the time contour in the path integral). 

2.2 The Full Lagrangian 

The interaction terms in the contour Lagrangian (2.2) involve quartic powers 
of x±- The momentum-space transformation (2.10) from x± to <Pr/a becomes 
in coordinate space 

X±(x) = J dy (p a (x + y , x)V a ±(y ), (2.12) 

where a is summed over R/A and V a ±(y ) is the Fourier transform of V a ±(k ). 
(See Appendix A.) The interaction terms (x±) 4 produce products of the field 
operators 4>r/a at four different times so that C c is not a true Lagrangian den- 
sity when expressed in terms of 4>r/a- To avoid products of non-simultaneous 
fields it is instructive to perform a heuristic Taylor series expansion. The idea 
is to express x± i n terms of the time derivatives of 4>r/a at the same x . A 
formal expansion of (p a (xo + yo,x) in (2-12) as a Taylor series in y gives: 

X±W (2.13) 

The symbol ~ is a warning that the formal series is meaningless but it serves 
to motivate what follows. The coefficients C e are defined as moments with 
respect to the time variable in (2.14) but they may also be expressed in terms 
of derivatives in frequency as in (2.15): 

CL= J dy ^fv a± (y ) (2.14) 

—00 
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r e _ HY d*V a± (k 
a± t\ dki 



(2.15) 



The explicit forms of V a ±(ko) shown in Appendix A contain Bose-Einstein 
functions of ko and are not infinitely differentiable at ko — 0. Since the actual 
value of these coefficients will not matter until the very end of the calcula- 
tions in Sec 6, there is no problem in regularizing the C e by introducing a 
fictitious chemical potential as discussed in Appendix A. At the end of the 
calculation the regularization may be removed. In the case of fermion fields 
the corresponding Vg±(/co) contain Fermi-Dirac functions and are not singular 
at k = 0. 

Using the series (2.13) makes all the terms in (x±) 4 have the same time. 
However the presence of an infinite number of time derivatives is difficult, and 
probably impossible, to quantize. The procedure instead will be to truncate 
the derivative expansion at a large order N, construct the Hamiltonian for the 
truncated problem, and later take the limit N— > oo. To do this, define the sum 
of the first N-l time derivatives as 

N-i u 



£=0 



so that if N— > oo then one expects <p± to approach x±- The truncation gives a 
Lagrangian density with quartic interactions (0+) 4 and (0-) 4 containing time 
derivatives up to order N-l. The method for constructing a Hamiltonian with 
higher time derivatives will be explained below. The essential feature is that 
there will be N-l canonical momenta. To construct the Hamiltonian requires 
inverting the relation between the highest time derivative and the highest 
momentum. This inversion is not manageable since the time derivative of 
order N-l is buried in the series for <f>±. There is a simple solution: Add an 
auxiliary term containing a higher time derivative of order N. The specific 
Lagrangian density to be used is 

t = (^ fi )(y^)-mU - h ^ QN J N A -A(0 + ) 4 + A(0_) 4 . (2.17) 



Now the highest time derivative is N'th order. Inverting the relation between 
the N'th time derivative and the N'th canonical momenta will be trivial. In 
order that h has no physical effects it will have to approach zero as N— > oo. 
The limiting behavior of h that will be used is 

h = hoN~ N . (2.18) 
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More general possibilities are discussed following (5.23). The field equation 
that follows from (2.17) is 

(d-d + m 2 )<fe = -4A £(-1)' C^al^W- ~ (2.19) 

where A — R and R — A. The highest time derivative in the field equation is 
2N; the highest time derivative of 4>r/a in the summation over £ is order 2N-1. 



3 The Full Hamiltonian 



With a finite number of time derivatives in C it is straightforward to construct 
the Hamiltonian for this system. The necessary generalization was developed 
by Ostrogradski in 1850 and is discussed in the mechanics book by E.T. Whit- 
taker [16]. For the Lagrangian density (2.17) there are canonical momenta 
, n^, . . . , ITj and canonical coordinates . . . , $^ where a = R/A. 

The capital letters II and $ will denote Heisenberg operators. (In Sec 4 and 
subsequently lower case 7r and will be used for interaction picture operators.) 
The canonical momenta are 

N—m q£ q r 

Momenta n*, Ii 2 a1 . . . , II^ -1 all involve the quartic interaction. Writing these 
out explicitly is tedious and not very useful. Only the last momenta has a 
simple form: IT^ = —hd N <j)a/dxQ. The canonical coordinates are the time 
derivatives of the fields: 

*i = ?w (i=l,2,...,N). (3.2) 



The Hamiltonian density is 



n= E (J« + + . . . + n?" 1 *? + n? f) v 

a=R/A OX 

-C(^,...,^,Q) (3.3) 

Note that this is only a Legendre transform with respect to the highest time 
derivative. This highest derivative is eliminated using d N (p a /dxQ = —U^/h. 
The explicit Hamiltonian density is then Ti = H + Hi where 
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^o= E (« + ^ + ... + n^)-JnX 

+ (V<^) • (V$\) + m 2 $^$\ - (3.4) 

ft 7 = A($+) 4 -A($_) 4 (3.5) 
It is very important that $± are functions only of the canonical coordinates: 

The nineteenth century development of higher derivative theories was clas- 
sical, but if the classical Poisson brackets are replaced by commutators the 
textbook treatment [16] becomes quantum mechanical. The canonical com- 
mutation relations are 

[Il™(t,x),^(t,$)} = U im 5 a p5\x-y). (3.7) 

All other equal-time commutators among IPs and $'s vanish. The Hamiltonian 
H — J d 3 xH generates a complete set of Hamiltonian equations for the 2N 
coordinates and 2N momenta: 

ifr t[H ^ ] &f =M (3 - 8) 

It is easy to check that this procedure really works. For each a there are 2N 
Hamiltonian equations (3.8). Of these, 2N-1 equations express the n and $ as 
time derivatives of . The last equation II 1 = i[H, U 1 } reproduces the field 
equation (2.19). 



4 Prelude: Quantization of C 



oo 



The construction of field operators that satisfy the Heisenberg equations (3.8) 
will be done perturbatively using the interaction picture with free dynamics 
specified by H . The corresponding free Lagrangian density is 

Co = {dM{d^ A )-m 2 <t> R <t> A -h ° d ™ ° d ^ A . (4.1) 



The presence of the order N time derivatives will be crucial in the next sections 
because it enables lower time derivatives of the fields to commute at equal 



10 



times. The ultimate strategy will be to compute matrix elements at fixed N 
and then let N — > oo. To understand why the correct Feynman diagrams 
emerge in this limit it is useful to first quantize the much simpler Lagrangian 
density 



A)0 = (<9 M R )((9>A)-m 2 0i?0A- 



(4.2) 



The canonical momenta for this are n R = d<j)A/dt and tia 
Hamiltonian density is 



d</> R /dt. The 



^oo = ^a-kr + (V0 fl ) • (V0a) + m 2 (j) R (j)A- 



(4.3) 



The equations of motion that follow from (4.2) are {d-d + m 2 )(j)R{x) = and 
{d-d + m 2 )(f)A{x) = 0. The solutions to the free field equations can be expanded 
in terms of spatial Fourier transforms. For later purposes this will be done in 
a box of volume V with periodic boundary conditions that make the wave 
vectors k discrete. The free particle energy is E — {m 2 + k 2 ) 1 / 2 : 



<f>A(t,X) 



(f> R {t,x) 



A ik-x—iEt I A —ik-x+iEt\ 



(&5 ^ik-x—iEt _|_ ^—ik-x+iEt^ 



Quantization is imposed by the canonical commutation relations 

1 



^r{x)Ar{v) 



%o=yo 



n A {x) : <p A {y) 



xo=yo i 



5 3 {x-x') 



(4.5) 



The equal-time commutators of [ttr, it a] and [<J)r, <Pa] are required to vanish. 
These can be satisfied by mode operators that obey 



b i>4\ = [ b Z>4\= 5 w/ 2E - 



(4.6) 



All other commutators vanish: \b?, cf, 



• ■ • • ••„••„, [b$,4\ = = 0. This im- 

plies the rather unusual result {4>r{x), 4>R{y)} = a{%) ■> 4> a{v)] — for any 
separation x — y. The only non-vanishing commutator of the two fields is 



[<Pa{x),< 



'R 



I 



d 3 k 



2E{2tt) 3 



e ik-(x-y) ^ e -iE(x -y ) _ e iE(x -y )^ {4.7) 



which will be important later. 
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Retarded and Advanced Vacuum States: Since the operators b A and b R com- 
mute, it would be conventional to define the vacuum state to be annihilated 
by these two. However it is more useful to define a vacuum state \0 R ) that is 
annihilated by the commuting operators b A and c A : 

g|0 fi )=0 if qe{b A ,c A }. (4.8) 



The state \0 R ) will be called the retarded vacuum. This state has an infi- 
nite norm and cannot be used to compute expectation values of the form 
(Or I O | Or). Instead it is necessary to introduce an advanced vacuum \0 A )- The 
advanced vacuum will be in the bra vector form (0a| to construct matrix ele- 
ments (Oa|C|0.r) and therefore it is natural to express its definition in terms 
of the operators that annihilate the bra state: 

(0 A |g = if qe {b R ,c R }. (4.9) 



The inner product of the two states can be chosen to be (CU| Or) = 1. 

Free Propagators for 4>r/a : The free propagators are vacuum matrix elements 
of time-ordered products. Since there are two fields 4>r/a the propagator is a 
2x2 matrix: 



A a/3 (x - y) = -i(0 A \ 



T{(p R (x)(f) R (y)} T{</> R (x)<f> A {y)} 



K T{(f)A(x)(f) R (y)} T{(/) A (x)(j)A(y)} 



\0r). (4-10) 



All four operator products are time-ordered, just as at zero temperature. The 
field operators and the states are all temperature- independent. The evaluation 
of (4.10) is easy. The relation (0^10^ = implies that (i) A RR = 0, (ii) A RA 
vanishes if x > y , and (iii) A AR vanishes if x < y . The relation 4>a\^r) = 
also implies (ii) and (iii) in addition to (iv) = 0. The no n- vanishing 

entries are thus 



A A r(x -y) = -i9(x - y o ){0A\<f>A{x)(f> R (y)\0 R ) 
Ara(x -y) = -i0(y o - x o ){0 A \<f>A(y)<f>R(x)\0 R ). 



(4.11) 



The matrix elements may be written as commutators since (0a\4>a(Pr\0 r ) = 
{®a\[<I>a,<I>r]\Qr)- The value of the commutator is given in (4.7). The Fourier 
representation for the propagator matrix is 



A af3 (x -y) = J 



d A k 



-ik-(x-y) 



A RA (k) 
A AR (k) 



(4.12) 
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where 

A AR (k) = l/{k 2 - m 2 + itk ) A RA (k) = l/(k 2 -m 2 - iek ) (4.13) 

are the usual free retarded and advanced propagators. This is precisely the 
same as (2.5) and will be the starting point of the perturbation theory devel- 
oped in the next section. 

5 Interaction Picture: Quantization of £ 

5.1 Field Operators in the Interaction Picture 



The full Hamiltonian (3.4)-(3.5) is expressed in terms of operators II and $ 
in the Heisenberg picture. Practical calculations in perturbation theory will 
require interaction picture operators it and 0, whose time evolution is given 
by Hq. The corresponding Lagrangian density is 

A) = (d l ,(f) R )(d >1 (f) A )-m 2 (f) R (f) A - h ■ (5-1) 

Quantization requires the canonical coordinates 

*«=-Qj^r (m = l,2,...,iV) (5.2) 



and the canonical momenta 



N-e 



where a = R/A. In terms of time derivatives, the momenta are 



< = (-l) N+1 ~ e h (t = 2, 3, . . . , N). (5.3) 

As usual the momenta all commute at equal times, the coordinates commute 
at equal times, and the only non-vanishing equal-time commutator is 

[ni(t,x),^(t,y\ = ±6*™6 3 (x-y). (5.4) 
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This implies that the first 2N-1 time derivatives of 4>a commute at equal times, 
the same for <p R , and that commutators involving time derivatives of <f>A an d 
of must satisfy 



d e <j) A (t,x) d m Mt,y) 



ro 



ifi+m ^ 2JV-l,4JV-3, 



i(-iy +N 5 3 (x -y) ite+m = 2N-\. 



(5.5) 



The cases in which i + m = AN — 3, 6N — 5, . . . are omitted because they are 
not independent of the case £ + m = 2N—1. They are related by field equation 
that follows from (5.1): 

= ( d . d + m 2 )^ + {-lfh°-^. (5.6) 



The field equation is solved by an exponential exp(?A; • x — iuijt) if ujj satisfies 

= (ujjf ~E 2 -h {u^ (5.7) 



where E = (m 2 + k 2 ) 1 / 2 . This equation has 2N roots. For every root ujj there 
is a root — O0j. For h very small, one root is real and close to E. This real 
root will be called u\. Thus — E + O(h). Since (5.7) is even, —lo\ is also a 
root. For h small, the other 2N — 2 roots are very large: \coj\ ~ /i -1 /( 2JV-2 ). For 
ft, ps jV - ^ this means that ss iV 1//2 for large iV. It is convenient to label 
the roots so that 002,003, . . . ,oon are all in the lower half plane. The roots in 
the upper half-plane are the negatives of these. The general solution to the 
field equation may be written 

I N 

4> A {x) = ^=EE ( b f% eik ' 3 ~ i0Jjt + c fk e- ik - s+i ^). (5.8) 



This form will be used for computing commutators but it disguises the self- 
adjoint property of <\>a- Since uj\ is real, = bf^. The other 2N-2 complex 
roots of (5.7) come in complex conjugate pairs. Thus for each one of the roots 
oo 2 ,oo 3 , . . . ,uo N in the lower half-plane the complex conjugate appears on the 
list — uo 2 , —003, . . . , —oon. For example, 002 is in the lower half-plane and suppose 
that a>2 = —^8- Then b^ = . This allows 4>a to be self-adjoint even though 
(5.8) conceals it. For the retarded field the solution is of the same form: 

1 N 

4 R (x) = ^=EE (b% e ik - x -^ + c% e - ik -^). (5.9) 
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The equal-time commutators (5.5) can be satisfied as follows. All the advanced 
operators, b and c A , must commute among themselves. All the retarded op- 
erators, b R and c R , must commute among themselves. The only no n- vanishing 
commutators are [bf, c R ] and [cf, b R ] for i—j. The value of this commutator 
will depend on the roots of (5.7). The only non- vanishing commutator among 
the mode operators is 

H><& = [%><& = s vQi s s,*> ( 5 - 10 ) 

where Qj is defined as the residue of the root ujj in the lower half-plane: 

(5.11) 



UJ 



2 - E 2 - hoo 2N 



The Qj satisfy two important identities 



N 

= Y,( UJ j) S Qj (sodd,s^2iV-l,4iV-3,6iV-5,...) (5.12) 

N 

1 = -2h^f N - l Q 3 , (5.13) 

3=1 

which are proven in Appendix B. These identities implement the commutation 
relations (5.5) as follows. Using the mode expansions and the commutators 
(5.10) gives 



ik-(x-y) N 

k v J=l 



If £ and m are both even or both odd, this vanishes trivially. If one is even and 
the other is odd, then £ + m is odd and the sum vanishes because of (5.12) 
except when I + m = 2N — 1, 4N — 3, . . .. For the case I + m = 2N — 1 the 
sum is given by (5.13) and has the correct value to agree with the canonical 
requirement (5.5). 



5.2 Free Propagators for (f) R / A 



In order for the usual time-ordered perturbation theory to succeed the free 
propagator must be the vacuum matrix element of a time-ordered product. 
As in Sec 4 vacuum matrix elements will be defined in terms of an advanced 
bra vector (0^1 and a retarded ket vector |0^). The natural manner of defining 
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these is through the mode operators b R / A and c R / A . For a fixed mode operator 
q, either q\0n) = or (0 A \q = 0. This partitions the mode operators into 
two disjoint sets. All the operators in each set are mutually commuting. The 
precise makeup of these sets is as follows: 



q\0 R ) = ifqe{c A ,b A X /R ,bf R ,...,bT} 

(0 A \ q=0 ifqe {cf , 6f , c A/R , cf /R , 4 /R } (5-14) 

The roles of b x and of Ci are the same as in Sec 4. The remaining operators 
are conventional in that for j > 2 the bj are destruction operators and the 
Cj are creation operators. Matrix elements of the form (0A\qq'\^R) vanish if q 
and q' are both 6's or both c's. If one is b and one is c the matrix element also 
vanishes if both are retarded or both are advanced. The only non-vanishing 
vacuum matrix elements of two mode operators are the following: 

(0A\b^\0 R ) = {Q A \b%c%\U R ) = 6ij Qj 6& for j = 2, 3, ... N. (5.15) 

As a result (OaIM^MWr) = and (0a\Mx)MWr) = 0. The 2 x 2 
propagator matrix is 

G af3 (x) = -i{0 A \T (Mx)MO)) |0 fl ), (5.16) 

and has vanishing diagonal components: 

I G RA (x)\ 
G aP {x) =\ RAK ' . (5.17) 



Gar(x) J 



The nonvanishing propagators are Gra(x) = — i(0A\T((f>R(x)(j>A(0))\0R) and 
Gar(x) = —i(0 A \T((j)A(x)(j)R(0))\0R). By translation invariance Gra(x) = 
Gar(-x). The result of evaluating Gar(x) using (5.15) is 



73 1 

G AR (x) = -i J e **Q x 9{t) (e-^-e^) 



r d 3 k - N 

~ % \ t^^^Qj «(«)^' + «(-t) e ^ . (5.18) 

J i^ 71 ) j=2 

The time dependence on the second line can be written exp(—iidj\t\). Since 
the ojj for j = 2,3, ... ,N have negative imaginary parts the second falls 
exponentially as t — > ±00. 
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Large N Limit: As A — > oo the first residue Q 1 -> 1/222. For j = 2, 3, . . . A 
the behavior is Q 6 -> -l/(2wj2V) ~ A~ 3 / 2 since u 6 ~ A 1 / 2 . Thus at large A 
the Gar propagator becomes the usual retarded free propagator (4.13): 

jim G AR {x) = -i6(t) J ^-y 3 e** (e~ iEt - e iEt ) = A AR (x) (5.19) 
and similarly for Gra{x) = Gar(-x). 

Time Derivatives of the Propagator: One of the unfamiliar aspects of the the- 
ory is that most time derivatives of the fields commute at equal times as 
stated in (5.5). These vanishing commutators are automatically built into the 
expression (5.18). The first time derivative dGAR(x)/dt contains no term pro- 
portional to S(t). However the second time derivative contains the term 

This automatically vanishes because of the residue identity (5.12). One can 
continue this process and show that d e GAR{x)/dt e contains no S(t) term for 
£ < 2 A. However for £ = 2 A there is a S(t) term whose coefficient is given by 
the inhomogeneous sum rule (5.13). 



5.3 Free Propagators for <p± 

A crucial test of the quantization is to calculate the temperature-dependent 
propagators 

G ab (x-y) = -i(0 A \T(Mx)Mv))\0R) M = ±) (5.20) 

since these are the naturally occurring quantities in perturbative expansions. 
The (ft± are defined in (3.6) as a sum of the first A — 1 time derivatives of (f>R 
and <pA- Thus 

JV-l Ql Qin 

G ab (x-y)= 7T^G^(x - 2/0) f-f)^ a C;, 

All time derivatives up to order A — 1 have been taken outside the time- 
ordering because all the equal time commutators vanish as discussed earlier. 
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Using the definition (2.13) of the temporal moments C l and C m gives 

G ab (x - y) = J dx' Q dy' Q X af3 V aa {x' Q )VM) (5.21) 



JV-l 



\x' Y d e \ f(y' o r d 



e,m=o 



t\ dx J \ m\ dyfi 



In the function X a p the yo derivatives can be converted to xo derivatives after 
which one of the two summations can be performed to give 

Xap = 2 £ - Vo, x - y). 

In this form X a p is recognizable as a Taylor series for G a /3 with the remainder 
term omitted. Since all the time derivatives of G a p up through order 2N — 1 
are continuous (with no 5(x — yo)), Taylor's theorem gives a specific form for 
the remainder R a/3 : 

X a p = G a(3 (x + x' Q - y - y' , x - y) - R a(3 (5.22) 

iA-y'o)^- 1 d 2N -'G a ^,x-y) 



Raf3 — 



{2N-l)\ d£ 



27V- 1 



where £ is a time in the range x — y < £ < x — y + x' — y' . For any value of 
N there is a £ such that (5.22) is an exact equality. By using (5.18) it is easy 
to see that the remainder vanishes in the large N limit: 

RaP (2N-l)\ ~ j2N~Ty. ~ h(2N-l)\ °" (5 ' 23) 



In the last step the behavior h ~ N~ N has been used. (More general h of the 
form h ~ N~ cN for c in the open interval < c < 2 will also work in that 
both the remainder (5.23) and the second line of (5.18) will go to zero.) The 
asymptotic vanishing of the remainder implies that 



lim X a/3 = Jim G afS (x + x' - y - y' ,x - y) 

= A a/3 (xo + x' -y -y' 0: x-y). 

because of (5.19). When this is substituted into (5.22) the time-integration 
variable can be shifted to a more convenient form as follows: 



18 



lim G ab (x -y)= dx' dy' A al3 (xo + x' - y - y' , x - y)V aa {x' Q )Vp b {y' ) 

TV— *oo J 

= J dx' dy' A af3 (x' -y' ,x- y)V aa (x' - x )Vp b (y' - y ). 

In momentum space this reads 

lim G ab (k) = VT(-k )A a p(k)Vp h (k ). (5.24) 

N— >oo 



Comparison to (2.6) and (1.1) shows that the right hand side is indeed the 
correct free thermal propagator: 



e^-^ ko n(A R - A A ) nA R -(l + n)A A 



Free Propagators for Mixed Combinations: 

It is sometimes necessary to evaluate propagators that involve one fundamental 
field 4>r/a and one (f>± field. The mixed propagator 

G ab (x -y) = -i(0 A \T((f> a (x)<i> b (y))\0 R ) (5.26) 



is evaluated in the same way. The time derivatives can be taken outside the 
time-ordering. Using the definition of C^ b as temporal moments of the thermal 
matrices Vp b leads to a single-variable Taylor series with the result 

N^L Gab ^ X ~ ^ = I dz o^( x o ~ z o,x-y)Vp b (z - y ). 



In momentum space this implies 

^G ab (k) = A a p(k)Vp b (k). (5.27) 



6 Time-Dependent Perturbation Theory 



All the necessary interaction picture propagators have been calculated. It is 
now straightforward to apply conventional operator perturbation theory in 
terms of interaction picture fields. The resulting Green functions will have 
the property that when N — > oo they give the correct perturbative results. 
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As usual the Heisenberg operators <3> a are related to the interaction picture 
operators <p a by 

$ Q (*,£) =U(0,t)(/> a (t,Z)U(t,0), (6.1) 



where the evolution operator is 

U(t, t') = exp[iH t] exp[-iH{t - t')} exp[-iH t']. (6.2) 

This may be expanded in terms of operators in the interaction picture: 

oo (_A)n K _ ~ n 

U(t, t') = J2 / dhdt 2 . . . d^TlH^Hjih) . . . HjiQ] . (6.3) 

n J 



71=0 



To compute the vacuum matrix elements of time-ordered products of the in- 
teraction picture fields requires a version of Wick's theorem. 

Wick's Theorem: To deal with time-ordered products of interaction picture 
fields it is convenient to introduce the generating functional 

S[j] = (O A |Texp(-^rf 4 a;j Q (a;)0 Q ( a ;))|O R ). (6.4) 

Define B(x) = j a {x)4>a{x) for brevity. Because the commutator of B with 
itself is a c-number the time-ordered exponential operator may be related to 
an un-ordered exponential operator that is multiplied by a c-number: 

Texp(—iJ d 4 xB(x)) =exp ( — ij d A xB(x)) 

xexp(- \J d 4 xd 4 y e(x -y )[B(x),B(y)}). (6.5) 

In the plane wave expansions (5.8) and (5.9) of 4>a and 0r, each of the mode 
operators either annihilates the retarded vacuum ket, q\0n) = 0, or annihilates 
the advanced vacuum bra, (CU|g = 0, as displayed in (5.14). This property 
provides a unique decomposition B(x) = Bi(x) + B 2 (x) where B 1 (x)\0 R ) = 
and (0a\B2(x) = 0. The vacuum matrix element of (6.5) is therefore given by 
the Baker-Campbell-Hausdorff formula: 

S[j] =exp(-i J d 4 x([B 1 (x),B 2 (y)]+e(x -y )[B(x),B 

Using [Bi(x), Bi(y)] = [B 2 (x), B 2 (y)] = this reduces to 

S[j] = exp(- J cPxcPyeixo-yoftB^xlBtiy)]). 
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The commutator is [Bi(x), B 2 (y)] = (0 A \B(x)B(y)\0 R ) so that the exponent 
can be written as a vacuum matrix element of a time-ordered product 

S\j] =exp(-i J d 4 xd 4 y (0 A \T(B(x)B(y))\0 R )). 

Restoring the definition B(x) = j a (x)(p a (x) gives the final result in terms of 
the propagator (5.16): 

1 t 

S[j] = exp ( -- J d A xd A yj a (x)G aP (x - y)jp{y)). (6.6) 



This is the generating functional for Wick's theorem. 

Sample Commutation: A typical computation is the full propagator for the 
Heisenberg fields $ a , a = ±. This is a 2 x 2 matrix 

G' ah {x-y) = -i(0 A \T(* a (x)* b (yj)\0 R ). (6.7) 
It can be written in terms of interaction picture fields: 

(Oa|c/(oo, -ooJlOij) 

This can be expanded in powers of the coupling strength, A, as G' = G^> + 
G^ + G^ + . . .. The lowest order term, G^°\ is the free propagator computed 
in Sec 5.3. As shown in (5.25) g£ } goes to the usual free propagator matrix A a b 
in the limit limit N^oo. The first order correction to (6.8) is the connected 
part of 

G^(x-y) = j d'z(0 A \T(Mx)[-X<l>X(z) + X^iz^My))^}. (6.9) 

Wick's theorem (6.6) allows the vacuum matrix element of the time-ordered 
product of six fields to be written in terms of vacuum matrix elements of three 
pairs of fields. The term gives 

-12A J d 4 z {0 A \T(<l> a (x)<l> + (z))\0 R ) (0 A \T(<l> + (z)<i> + (z))\0 R ) 

x{o A \T(<i> + (z)My))\OR)- 

After a similar evaluation of the 0i term the result is 

G$(x -y) = I2i\ J d 4 z ( G a+ (x - z)G ++ {z - z)G +b {z - y) 
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- G a 4x - z)G__(z - z)G_ b (z - yj) (6.10) 

In momentum space this is 

G$(k) = 12t\ J -0- 4 ( Ga+ (k)G ++ (q)G +b (k) - G a -(k)G-4q)G- b (k)). 

Using the result (5.25) the large N limit is expressed in terms of the usual free 
propagator A ab 

Hm Gff(fc) = l^A J A a+ (k)A ++ (q)A +b (k) 

- A _(A;)A__(g)A_ 6 (A;)). (6.11) 

This is the correct result for the one-loop correction to the thermal propagator 
matrix and agrees with the standard answer [6-10]. From this example it is 
easy to see that one can compute order by order in the interaction picture 
and in the limit N — > 00 the result will coincide with the usual perturbative 
expansion. 



7 Discussion 

The purpose of constructing the Hamiltonian formulation of finite temperature 
field theory in terms of <pR and 0^ is not to obtain another derivation of the 
Feynman rules for the n-point functions. On the contrary, the Feynman rules 
were developed in order to show that the temperature-dependent Hamiltonian 
H has the same physical content, i.e. exactly the same perturbative expansion 
of the Green functions, as the usual formulations. 

Any quantum mechanical system allows both wave and particle descriptions. 
Finite-temperature field theory to date has been quite successful with the 
wave aspects in that the n-point functions can be calculated perturbatively. 
The path integral formulation places particularly strong emphasis on the wave 
properties since the Feynman rules may be obtained there without ever intro- 
ducing the notion of a particle. However the particle aspect of the theory can- 
not be put off indefinitely. The problem is that when particles are in a thermal 
background their self-energies become temperature-dependent. Consequently 
a propagator that normally has a pole at k$ = m 2 + k 2 will instead have a 
pole at some ko that is temperature-dependent. Such a pole corresponds to a 
thermal quasiparticle. Important evidence that this pole is the proper way to 
characterize the quasiparticle comes from the theorem of Kobes, Kunstatter, 



22 



and Rebhan [17] that in any gauge theory the location of this pole will be 
gauge invariant. 

The operator formulation presented above is a preliminary step in under- 
standing the quasiparticle poles. As discussed in the Introduction, the full 
finite-temperature retarded propagator D' R (k) has a quasiparticle pole in fc . 
As shown in Sec 6 this same propagator may be computed as 

D' R (x -y)= lim \(0 A \T(* A (x)* R (y))\0 R ). (7.1) 

Because <£> are Heisenberg fields with their time dependence given by the 
temperature-dependent H, $(t,x) = exp $(0, x) exp(— iHt) and conse- 
quently the Fourier transform of (7.1) is 

D' R (k ,x-y) = hm (O A \<S> A (0,x)- l—* a (0,jft\O R ). (7.2) 

N^oo fc Q — H + IT] 

A future publication will show how the poles in ko come from isolated eigen- 
values of H and that the branch cuts in ko come from continuum states of 
two or more quasiparticles. Thus H determines both the time evolution of the 
fields and the energy spectrum of the particles in direct analogy with zero 
temperature. It is surprisingly easy to adjust H$ so that its spectrum coin- 
cides with that of the full Hamiltonian H. The shifted Hq will describe free 
quasiparticles. This is parallel to the zero-temperature shift that make the free 
Hamiltonian contain the physical mass rather than the bare mass. At finite 
temperature this leads to a quasiparticle resummation program. Although the 
above discussion only treats spinless fields, a future publication will extend the 
method to the more physically interesting case of fermions and gauge bosons. 
The ultimate purpose is to obtain probabilities for quasiparticle processes that 
are gauge-invariant and infrared finite. 

The methods used in Sec 3 may be useful in other problems. It often hap- 
pens that effective actions contains products of fields a different times. The 
Braaten-Pisarski effective action [18] is such an example. Sec 3 and the later 
development shows how to expand the action in a finite number N of time 
derivatives, apply canonical methods to construct the Hamiltonian or the 
energy- momentum tensor, more generally, and quantize the result. 
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A Diagonalization of the Propagator 



The matrix V a b(k ) that diagonalizes the free thermal propagator in (1.1) was 
shown by van Eijck, Kobes, and van Weert [13] to be 



f V R+ (ko) V R _(k )\ _ I e^n(k )/a(-k ) e° k °n(k )/a(-k ) 
K V A+ (k ) V A -(fo)j V a ^°) e^a(k ) 

where a(k ) is arbitrary. The matrix U a b(k ) that occurs in (2.7) is 



,(A.l) 



( U R+ {k ) U R _(k ) 
K U A+ (k ) U A -(ko) 



a(ko) -e° k °a(k ) 
e^n(k )/a(-k ) -e° k »n{k )/a{-k ) 



(A.2) 



The expansion coefficients C l ab are defined in (2.15) as the £'th derivatives of 
Vab(ko) evaluated at k = 0. The Bose-Einstein function in (A.l) does not 
have finite derivatives at ko = 0. If a(ko) is chosen to have a pole at ko = 
then the first row of (A.l) will have finite derivatives but not the second row. 
It is easy to evade this problem as follows. As shown in Sec 6 it is only at the 
end of the calculation, after calculating Green functions and taking the limit 
N — > oo that V ab (k ) has to be (A.l). In the intermediate stages V ab (k ) can 
differ from (A.l). One simple possibility is to introduce a fictitious chemical 
potential fi and set n(ko) = l/[exp((3(ko — //)) — 1]. In intermediate stages 
of the calculations the derivatives of n(ko) will give inverse powers of /i. At 
the end of the calculation, the limit N — > oo sums all the inverse powers of 
fi to give back l/[exp((3(k — //)) — 1]. Thus one obtains a perturbative series 
containing l/[exp((3(k — //)) — 1] in the propagators. At that point it is safe 
to set /i — > to obtain the true theory. Note that since the particles under 
consideration are massive there are no infrared divergence. 



B Residue Identities 



In the interaction picture the differential equation (5.6) for the field operators 
leads to the frequency equation = u 2 — E 2 — hu 2N . For h very small there are 
two small roots of the form cu ?a ±(E — hE 2N ~ l /2). There are N-2 large roots 
uj w ^-V(2JV-2) exp ^ 7r £/(iv - 1)], where i = 1, 2, 3, . . . , 2N - 2. The particular 
root with £ = N — 1 is negative real; that with t = 2N — 2 is positive real. 
Thus of the total 2N roots, there are two positive real, two negative real, N-2 
complex roots in the lower half-plane, and N-2 complex roots in the upper 
half-plane. If the equation is modified to = uj 2 + irj — E 2 — h(uj 2 + ir/) N for 
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T) infinitessimal then the two positive roots acquire small negative imaginary 
parts and the two negative roots acquire small positive imaginary parts. Then 
there are N roots in the lower half-plane and N in the upper half-plane. 

For s = 1, 3, 5, ... , 2N — 3 the following integral is convergent but vanishes 
since the integrand is odd: 



oo 

/ 



duj 



oj 2 +ir ] -E 2 -h( y u 2 +ir ] ) N ' 



The contour can be closed in the lower half-plane. There is no contribution 
from the semicircular contour at infinity since s < 2N — 3. There are N roots 
ujj in the lower half-plane. By Cauchy's theorem and the definition (5.11) of 
the residues Qj the integral gives 

N 

= J2M S Qj (s = 1, 3, 5, . . . , 27V - 3). (B.l) 

.7 = 1 



This result can be extended to higher values of s as follows. Keep s odd but 
in the range 1 < s < 2N — 5. Then multiply each (ujj) s Qj in (B.l) by 

E 2 = (a;,) 2 - h{u 3 ) 2N . (B.2) 

The sum with (ujj) s+2 vanishes by the previous result. Hence the sum with 
(cUj) s+2N also vanishes, where 2N + 1 < s + 2N < AN — 5. No conclusion 
has been drawn about the sum with (ujj) 21 ^^ 1 . The same technique can be 
repeated. The complete result is 

N 

= Y^M'Qi (sodd,s^2iV-l,4iV-3,6iV-5,...) (B.3) 

To evaluate the sum at the first exceptional value, s = 2N — 1, use the identity 



U 2 _ E 2 _ hu 2N ~ E Qi U 2 > (B-4) 

which defines the residues Qj. Set oo = and multiply both sides by E 2 : 



N O 

^3 77.2 



3=1 ^ 
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For each term in the sum use (B.2) and the previous result (B.l) to get 



N 

-l = 2hJ2QA^) 2N - 1 (B.5) 

3=1 

To evaluate the sum at the next exceptional value of s, multiply both sides of 
(B.4) by E 4 , apply d/du 2 , and then set uo = 0: 

For each term in the sum use (B.2) and the previous results (B.l) and (B.5) 
to get 

N 

-l = 2h 2 Y / QM) 4N - 3 . (B.6) 

3=1 

The next case is 

N 

-l = 2h 3 Y,QMf N ~ 5 - (B.7) 

3=1 
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